Scattering amplitudes in non-Fermi liquid systems 



A. Koutouza^, F. Lesage^ and H. Saleur^^'^'* 
^Department of Physics, University of Southern California, Los Angeles, CA 90089-0484 
^ Centre de recherches mathematiques, Universite de Montreal, 
CP. 6128, succ. centre-ville, Montreal, P.Q., Canada, H3C-3J7 
Physikalisches Institut, Albert-Ludwigs-Universitat, D-79104 Freiburg, Germany 
... ^ Service de Physique Theorique, CEN Saclay, Cif Sur Yvette F91191 

■ (Dated: February 2, 2008) 

o ■ 

' By a mix of form-factors and analyticity techniques, we determine some fundamental scattering 

^S) , amplitudes in non-Fermi liquid systems. These include the reflection and transmission amplitudes 

, ^ . for Laughlin quasiparticles at a point contact between two g = \ Luttinger liquids, and the e — > e, 

^ l ' Iqp — > Iqp and e — > Cedge at a point contact between a Fermi liquid and a (j = i Luttinger liquid 

, (all liquids chiral or not) . These results are obtained in closed form, and give rise to rather simple 

■ expressions for the probabilities of the most basic processes of non Fermi liquid physics at these 
04 ' special values of the couplings. 
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PACS numbers: 11.15.-q, 73.22.-f 



I. INTRODUCTION 



■ The purpose of this paper is to address the calculation of scattering amplitudes in systems exhibiting non-Fermi 
\ liquid physics. In such systems, scattering processes involving electrons (more precisely, the quasiparticles adiabatically 
connected with the bare electrons in the Landau Fermi liquid theory, with charge e and spin 1/2) can give rise to the 
appearance of other kinds of quasiparticles. A characteristic example of this phenomenon (to be discussed in details 
Q . below) occurs when electrons are injected from a Fermi liquid into a Luttinger liquid, and can give rise to various 



combinations of Laughlin quasiparticles and quasiholes I*!. Another example occurs in the spin 1/2 two channel 
Kondo problem: at the low energy fixed point, the amplitude for an electron to be scattered into an electron (or any 
finite combination of electrons and holes) is actually zero, hence giving rise to ' unitarity puzzles', which are solved 
by recognizing that the scattering takes place entirely into the spinors sector Pj . 

■ The determination of these amplitudes is conceptually a very important question, as one may argue they are the 
' essence of non Fermi liquids physics. Technically however, it is a very difficult one, which has been tackled mostly 

■ in perturbation theory around the high energy fixed point 0. The purpose of this paper is to present some non 
"nI ■ perturbative attempts at calculating these amplitudes in a case where the problem is integrable - and especially 

simple. 

Our initial attempts were based on the technique of form-factors and this paper builds on earlier works in this 
area As we will see, the answer to the physical questions of interest involves the determination of correlation 

"j^ functions in a massless field theory with a boundary interaction. In the paper Q it was shown that such correlation 
functions could be determined with remarkable accuracy in the case of operators with no anomalous dimension, like 
I . the current and the stress energy tensor. In the paper 0, preliminary attempts were made to determine similar 
'■^ ■ correlators in the case of vertex operators with non trivial, anomalous dimensions. From a technical point of view, 
Ch , the present work is a follow up of that paper, and will use similar regularization techniques. It will turn out however 
O ■ that for the specific questions asked - for instance, what is the probability that an electron from a (chiral) Fermi 
' liquid tunnels into a (chiral) Luttinger v = ^ liquid as an edge electron, the form-factor results do not, in practice, 
K*" ■ give reliable results for all scales of energy. We will thus have to make considerable use of another non peturbative 
("analytic") approach pioneered by Chatterjee and Zamolodchikov 8] to obtain what will turn out to be results in 
. closed form. 

' Two main applications are presented: the problem of an impurity in a one dimensional Luttinger liquid, and the 
problem of a point contact between a Fermi liquid and a Luttinger liquid. These two problems are related with 
the boundary sine-Gordon model, with which we assume the reader to be familiar. Mostly for technical reasons, we 
restrict in this paper to the case g = \- Although this is a "free fermion theory" , the amplitudes of interest involve 
operators which are not local in terms of the fermions, and whose correlators are already quite complex. The physical 
features are not expected to depend strongly on g, so we expect our results to shed light on the general situation; 
moreover, g = ^ is of direct relevance to the case of a contact between a Fermi liquid and a. v = edge :9;3. 

In the case of the impurity in a Luttinger liquid, or, equivalently, tunneling between two identical chiral Luttinger 
liquids (section 2), we discuss both the reflection and transmission amplitude for a single Laughlin quasiparticle. In 
the case of tunneling between a chiral Luttinger liquid and a Fermi liquid (section 3), we discuss the e — > e, the 
Iqp — *■ Ipq and the e — *■ Cedge amplitudes. Here Cedge denotes the edge electron, which can also be considered as the 
bound state of three Laughlin quasiparticles. 

Exact expressions are obtained in all cases. For the ease of the reader, we gather these expressions here. In the 
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FIG. 1: Tunneling between two identical chiral Luttinger liquids. What we call reflection corresponds to a Iqp remaining 
within the same liquid, and what we call transmission amplitude corresponds to a Iqp jumping from one liquid to the other. 
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Some technical issues are relegated to the appendices. 

To conclude this introduction, we would like to point out that the mere questions we want to answer have to be 
defined very carefully, a task we do not completely undertake in this paper. Indeed, in the massless theories we are 
discussing, there are ambiguities in the set of scattering states one uses as a "physical basis" . For the one dimensional 
Fermi liquid, the choice of electrons and holes is canonical, but for general Luttinger liquids, this is not so. For instance, 
for the g = \ case, a choice based on charges i (Laughlin quasihole) and — 1 (edge electron) quantas has been studied 
in great details Other choices, eg using charges i and — ^ (Laughlin quasiparticle) quantas, should be possible 
as well: none appears more fundamental than the other, but the choice has to be carefully specified whenever one 
discusses, for instance, unitarity. Whatever the choice, the vertex operators (in the bosonized formalism) one uses 
as "creation/annihilation" operators of these quantas are not one particle operators in the usual sense. They do not 
simply add or subtract a single quasi-particle from a many particle state, but have a more complicated action due to 
their non trivial commutators. This can make the mere definition of a scattering amplitude quite non trivial |lOl |. 



II. IMPURITY WITHIN A LUTTINGER LIQUID 



A. Generalities 



We consider a problem with two identical chiral Luttinger liquids and tunneling between them induced by a gate 
voltage - the set up is represented on the figure ^ 
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The Lagrangian reads 

L = -^dxcjyiidt - da;)<j)i + -^d.x(f)2{dt - d.^)(f)2 + v5{x) ^^[^/^{(tyi - ^2)], (3) 
47r 47r 

where the fields 4>i, (f)2 are associated with each of the two Luttinger Uquids. By forming hnear combinations 

01+02 , 01 - 02 

the field 0' decouples, while the dynamics of the field is determined by the Lagrangian 

L = -^dx(t>{dt — dx)(f) + vS{x) cos(\/2i^0). (5) 

The coupling constant v describes the strenght of the tunneling process. The normalization is such that one has the 
propagator (T the time ordering operator) 

(T0(a;, t)0(O, 0)) = -^^ ln(a; - t + iesign{t)). (6) 

It is then convenient to fold this problem onto the half line x E [—00, 0], with action (here $ = (f>{x, t) + cf>{—x, t) is a 
non chiral boson) 

^=^y° dx[{d^^f)-{dt^f]+vcos[^^^m. (7) 

The boundary interaction has physical dimension d = v. In the following we consider mostly the case v = ^. 

The two amplitudes we wish to consider are what we call the reflection and transmission amplitudes for Laughlin 
quasiparticles. The creation operator of a Laughlin quasiparticle in the first Luttinger liquid is proportional to e"'^^'^^ 
so the reflection amplitude (corresponding to a Laughlin quasiparticle going through within the same chiral Luttinger 
liquid) is 

_ ^rp giv^0i(a:2>O,t2)g-i\/i^</>l(a:i<O,ti)j 

(T (,i\/W24>{x2>0,t2) ^-iy/U/24'{xi<0,ti)'^ ^g-j 



{x2 ~ xi — t2 + ti + ie sign{t2 — ti))'^/^ 

while the transmission amplitude (corresponding to a Laughlin quasiparticle jumping from one chiral Luttinger liquid 
to the other ) reads 

_ frjn gi\/!702(a;2>O,t2)g-i\/^</>i(a:i<O,ti)^ 

_ i (J' g-i-/;7720(2;2>O,t2)g-i-y/lV20(a:i<O,ti)^ ^g^ 

{x2 — xi — t2 + ti + ie sign{t2 ~ ti))"/'^ 

The main technical difference between these two amplitudes is the fact that exponentials of the field have opposite 
signs in the refiection case and identical signs in the transmission case. Restricting to = ^, the two objects we have 
to determine are thus 

5fl^T(l,2) = (T exp(±i0(a;2 > 0, t2)/2) exp(-z0(a:i < 0,ti)/2)). (10) 

To proceed, we first use the mapping to the folded theory on the half line. A field at a; < is an 'in', right moving 
field, while a field at x > now becomes an 'out', left moving field. The new correlators are then given by 

5fl,T(l,2) - (T exp(±i0i(-X2,i2)/2)exp(-z0fl(xi,ti)/2)). (11) 

The general approach to compute such correlators is the massless form- factors technique 0, @ . Let us therefore 
recall some basic facts about intcgrable boundary field theories. The hamiltonian is diagonalized using a quasiparticle 
basis, with left and right massless particles e = ±p = e*^, carrying a label ± corresponding to soliton or antisoliton. 
The bulk scattering of these quasiparticles is S = —1. Their boundary scattering is 

g/3-/3f5 _ j 

where e^^ oc f^. 
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The form-factors of the vertex operators can in principle be obtained by taking the massless hmit of the formulas 
in Subtleties arise in taking this limit however. Let us first recall the massive result. Restricting to \a\ < 1, the 
basic formula in reads, setting Ua — exp(iQ;$), 

(0|[/a(x,t)|g2»...6'l) + - + -~-~ ^ (-_^„(„-l)/2 /^sinTTOf 



(0|C/a|0) ' ' \ i 



2n 



exp j — im^^(i cosh^fc — xsinh 0^) j exp j a^^(^^„+fe — 
\ fc=i / V fe=i 

ni<fc<,<nSinhg^sinh ^"^-/"+- 
ni<fc,,<„cosh^^ 



(13) 



The vacuum expectation value itself is obtained using techniques of and is of the form < 0|[/c(|0 >= (am)" C(a) 
where a is a UV cut-off and C(a) a numerical constant whose expression is exactly known. Notice that the form 
factors are neutral, that is they are non zero only acting on states whose total solitonic charge vanishes exactly. 

The massless limit is obtained by letting rapidities go to infinity setting 9 = zLOq ± (3 where 6*0 — > oo, and letting at 
the same time the mass go to zero with me^" — > 2 (the choice of this constant is non universal, and can be absorbed 
in a global shift of the rapidities). For example, let us take the massless limit of the form- factor by choosing p of the 
n particles with charge -|- to be R movers, and g to be i movers (with p + q = ji) and similarly, u particles among 
the n particles with charge — to be R movers, v particles to be L movers (with u + v = n). As 9o oo, the leading 
behaviour of the ratio ifT^ behaves as 

exp[a0o(p -q-u + v)] ^s^pv^eouq 

This can be rewritten as e^^^" -(a+q-v) ] taking the prefactor (0|?7q|0) into account, it thus follows that the form- 
factors all vanish in the massless limit. The calculation of correlators is thus bound to be a complex task: what 
happens is that, while the form-factors vanish in that limit, the integrals over rapidities all diverge at low energy 
due to the presence of infinitely many soft modes, an "IR catastrophe" , making the contributions to correlators of 
undertermined form x oo. The proper way to proceed would be to keep track of this divergence carefully, by making 
calculations at finite mass, and then considering the limit when the length scales get much smaller than the inverse of 
this mass. It is not clear however how many form-factors would have to be taken into account in order to get accurate 
results - probably an infinity - and this procedure clearly is not efficient. 

The procedure we will use below bypasses this difficulty, and consists in working directly within the massless theory, 
essentially by considering ratios of correlators, which gets rid of the overall normalization problem. It turns out that 
the ratios \i'6\ have different behaviours depending on a and the balance of charges. In particular, the limit of the 
ratio is always finite when the charges are individually balanced in the left and right sector. In that case, the ratio 
factorizes, and one can individually define L and R versions as (setting y^.'^e^P-^ — gM0-R,i.) 



(0|K.|0) ' ' \ I 

X exp (^-i{t T ^) ^ ^^""^ ^(/5«+fe - /5fc) j 

^ ni<fc<,-<„sinh^sinh ^"+--^"^' 
ni<.,,<„cosh^^±|^ 

For I a I < i, this is in fact the only case where the ratio has a non zero limit. Our procedure thus seems well defined 
when |a| < i. When i < |a| < 1, some of the ratios diverge in the massless limit, and handling out the form- 
factors series is thus probably even more difficult. We have not tackled this case, but see the appendix for analytical 
techniques. 

The case \a\ — ^ is special. In addition to the case where L and R sectors are separately neutral, the ratio H13() 
also has a finite limit when the balance of charges is +1 in the R sector and —1 in the L sector. For instance 



while the same expression for y_i/2 vanishes in the massless limit. 
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Finally, the case |a| = 1 is also special, even in the massive case. This is because sinaTr vanishes, while the constant 
C(a) diverges. The best way to understand what happens is to make connection with fields with a — I and fermion 
operators in a free fermion theory. In the massless limit, one finds that V^^'" oc i/jR^Lii^ssp. tp]^ ^). This is discussed 
in more details below. 

The normalization of multiparticle states in (|14|l is < f3\(3' >= 2ttS{(] — (]'). We denote below the creation operators 
Z\°'j^{f3) such that eg .Z|''^(/3)|0 >= |/3 ^. Left and right creation/annihilation operators anticommute. 

The boundary interaction can be conveniently handled by moving to an Euclidian description, such that the impurity 
interaction becomes a boundary interaction, the boundary lying along the imaginary time axis. One then makes a 
a Wick rotation, and quantizes the theory along the boundary axis instead. The coordinate along the boundary, y, 
is related with the time hy y = it. The whole effect of the boundary interaction is then taken into account by the 
boundary state whose expression in the integrable basis is extremely simple: 



\B >= exp 



Ztt 



Here, the K matrix is simply related with the R matrix and given by 



K++ = K— - 



e/3-/3B + 1' 



K^- =K-+ = 



eP-pB + 1 



(16) 



(17) 



B. The transmission amplitude 



A calculation involving up to two particles gives then the following expression, where we neglected some trivial 
normalization factors (note that the L particles are only annihilated by the L field and the R ones by the R field) 

5t(1,2) = (exp(-i(/)L(-a;2,y2)/2)exp(-i0K(xi,?;i)/2)) 



/dl3 
— exp [e'^{xi -X2+ i{yi - 2/2))] — 



df3i d(32 exp [(e^i + eP^){xi ^X2+ i{yi - 2/2))] 



2ti 2-k 
1 



1 



1 + e''! 1 + 



„/32-/3i 



cosh^ ^ 



1 + e^^i 1 + ef^^ 



(18) 



where /? — P — Pb, and the dots denote processes involving three and more particles in the expansion of the boundary 
state. As mentioned brielfy earlier, this integral is divergent at low energies /3i — > —00 due to the proliferation of low 
energy particles, and the expression as it stands cannot be used. Integrals converge at high energy since X2 > 0,xi < 0. 

To proceed, we use a technique experimented in this context in 0. Denote ps = e''^ and observe that the low 
energy limit is also the limit v — > oo(pB 00), where the field (j) sees Dirichlet boundary conditions, and thus we do 
know that the correlator should go as 



5t(1,2), 



oc 



{x2 - xi + i{y2 - yi)y^* 



corresponding to a totally diagonal scattering. The trick to obtain finite results will thus be to consider not gx itself 
but the ratio 5t(1, 2)/gT{l, 2)pa=oo. Writing to first non trivial order gri^, 2) = 1 + /i + . . ., we thus consider the 
new approximation to the correlator 



5t(1,2) 



I + hips) 



{x2 - xi+ i{y2 
and expand the ratio to first order (ie, two /3 integrals) 

1 



5t(1,2) 



{X2 - xi+ i{y2 - yi))^/"* 



2/i))1/4 l + /i(oo) + ... 
[1 + /i(pb)-/i(oo) + ...] 



(19) 



(20) 



This procedure happens to give finite integrals, as is easily seen at this order, since the subtraction of /i(oo) cancels 
the divergence of the first integral in 118|l . It follows that, to leading order, the full correlator of interest reads, after 
the Wick rotation 



Gt(1,2) 



1 



{x2 - xi+ i{y2 - yi))^/^ 



1 

1 



du 



Jo Pb +u 



.^-u(x2-xi+i(y2-vi)) _|_ 



(21) 
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After a few simple manipulations using the Laplace representation of the square root prefactor, we end up with 



Gt{p,Pb) 



1 

1 - - 

TT 



p/Pe 



du 



, Pb 
1 u 



-1/2 



(22) 



In the low energy (that is, low energy of incident particles) limit p/pB ~> 0, corresponding to total transmission, 
Gt oc p~^/^(1 + est—), as expected. In the high energy Hmit p/pb ^ oo, Gt oc p^/^(l + cst\n{pB /p))- This latter 

Pb 

result clearly does not make much sense, as one expects Gt to vanish in this limit. 

A higher order calculation is perfectly possible, and would extend the degree of validity of the form-factors result 
towards high energies, while giving an essentially exact result at low energies. For any finite order however, the 
form-factors result seem to always become unreliable at sufficiently high energy, and behave like the first order in that 
limit. Rather than pursuing this direction (which might be the only one avaialble for other values of g), it is better to 
observe that the quantity Gt(1, 2) can be evaluated in closed form. This follows from a chain or arguments relating 
it to a very elegant calculation of the one point function of the spin operator in the Ising model with a boundary 
magnetic field carried out by Chatterjee and Zamolodchikov Q. The same kind of argument was used in the study 
of the Friedel oscillations in T?]. There, it was shown that that 



ST = (exp 



{a{x))h{<7{x))c 



(23) 



where (cr)^ denotes the one point function of the spin operator in the Ising model with boundary field h oc v (the 
right hand side does not depend on y when the two vertex operators are inserted at the same (imaginary ) time) , and 
X is the distance from the boundary. This is similar but different from the usual result in the bulk [Tslflql 



|sin($(z,z)/2)sin($(0)/2))|0) 
I cos($(z, z)/2) cos($(0)/2))|0) . 



(24) 



This result immediately extends to the case where the right and the left field are inserted at different points, since 
the correlator depends on X2 — xi + i{y2 — j/i) only: 



gr = {exp{-i(f>L{-X2,y2)/2)eiipi~i(l)R{xi,yi)/2)) 

= {cr{{x2 -xi+ i{y2 - yi))/2)h{cr{{x2 - xi + i{y2 - yi))/2))c 



(25) 



Using form factors, one can expand the one point function of the Ising spin operator as 

r-oo " 



{<j{x)), 



E: 

n=0 



1 



n 



,f3B- 

tanh 

27r 2 



Jl ( tanh 



(26) 



and one can also check the factorization l|23|l of g directly using the form-factors expansion as described in [T^.The 



case of infinite magnetic field corresponds to /i — s- cx) and since pB 



inh it also corresponds to ^ oo in 



this formula. In that case the result is known to be from boundary conformal field theory. Recall on the other 

hand the result that 



{a{x))hOC x^'\r^ eP-- Mpbx) 
V 47r 



(27) 



one thus find that 



Gt(1,2) oc e^--'^ IMpbx/2), 

V 47r 



X2- xi+ i{y2 ~ yi) 



(28) 



where the factor 1/2 came up in the argument of the Bessel function since in the initial Ising case, our variable x 
would be twice the value in (|27|l . We can now use the integral representation 



Ko{z) 



V2 Jo Vi 



^1/2 



dt 



to write 



Gt(1,2)(x 



'XPj 



-1/2 1 



PbJ 



-1/2 



dp 



(29) 
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giving rise to the stunningly simple momentum dependent transmission amplitude 

-1/2 



(30) 



Finally, we perform a shift (in the rapidity variable) of the integration contour to get back to the theory in real time, 
quantized with x as the space axis, and the interaction at the origin, which amounts to p ~f ip: 



Gr{p)oc (i + ^ 

PB 



-1/2 



(31) 



Note that this result could also be obtained by quantizing the theory in the crossed channel. In that case, the 
boundary interaction appears because asymptotic states are of the form \/3)r + R{f3)\(3)L, and leads to expressions 
with the shifted contour, and K replaced by R. This correspondence is discussed in details in |^. 

At low energy, Gt oc p^^/^(l + cst^) and at high energies Gt oc P^^^'^ ^J^- this is the expected behaviour. 
Indeed, in the limit of small pB the perturbative expansion should be linear in w cx y/ps, and one thus expects the 
expansion to go as \/pbIp- On the other hand, at large pb-, one approaches the IR fixed point (total transmission), 

along the operators cos$ of dimension d — 2, and ^9$^ [l/3|(here, ^ — (pR — (pL is the dual of the field 0, and the 

theory is defined on the half line). The coupling of these operators must therefore have dimension —1, and thus be 
proportional to since v has dimension 1/2. The perturbative series starts linearly in the coupling, that is linear 
in p/pb, as we just found. 

It is interesting to recall that the Bessel function Kq only admits an asymptotic expansion in powers of the argument 



E 

.fc=0 



r(A: + l/2) 



(2z)fe fc!r(-fc + 1/2) 



In our problem, this corresponds to the expansion near the IR fixed point (low energy of incident particles, large pb), 
which is thus only asymptotic for the correlator (in fact, the expansion of Kq is good only up to exponential terms 
e~^. Here this corresponds to terms e^^/^^/" \ which are non perturbative in term of the IR coupling constant l/w^)- 
All these difficulties disappear however in the momentum dependent amplitude, whose expansion is now convergent, 
thanks to the Laplace transform. 

Let us get back finally to the expression for {a{x))h- The leading order in h is proportional to hx^^^ \n{xh^). In 
perturbation theory meanwhile, the leading order involves an integral over the boundary (in the Euclidian formalism) 



, ^ 3^ 2)1/2 which must be regularized by the introduction of an IR cut-off. The situation is similar to the massive 

bulk Ising model or the bulk sine-Gordon model at = in (see and referecnes therein). An additional renormal- 
ization is required, leading to the replacement of the cut-off by the inverse of the mass - here, by the inverse coupling 
constant 1/pb- Alternatively, one could instead consider the correlator at finite temperature. 1/T then acts as an IR 
cut-off, and the scaling limit of (a) is entirely well defined. Taking the limit T — > then gives the result mentioned 
above jT3 |. 



C. The reflection amplitude 



The reflection amplitude is an intriguing quantity. Let us think of it first by considering the vicinity of the IR fixed 
point, which is is approached along the operator cos$. In an IR perturbation theory, the average of the quantity 
exp [«(0i? — (/)l)/2] will vanish to every order, because a charge neutral exponential can never be formed. Therefore, 
if the reflection amplitude happpened to be non zero, if would be entirely due to non perturbative effects! In fact, we 
think this amplitude is exactly zero. This can be argued more solidly by considering the UV point of view. The form of 
the IR perturbation shows that, right at the UV fixed point, the field at the boundary obeys (t>R — (phi^) = (2n + l)7r, n 
an integer. It follows that the quantity exp [i{4>R — 4>l)/2] changes sign when one goes from one well of the boundary 
potential to the next. Since a proper conformal invariant boundary condition involves a sum over all wells, the average 
(exp [i{(f>R — 4'l)/2]) has to vanish right at the fixed point. This observation can be made more formal by appealing 
to the formalism of conformal boundary conditions (notations follow those of [l9|: a„ are the usual bosonic modes; 
\w, k > the zero modes). The Dirichlet state is given by consideration of the IR fixed point (radius r = 



1 



V2t 



E 

k— — oc 



E 



|0,fc) 



(32) 
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The Neumann state compatible with this is 



^-2inrw<S>o g^p 



.n=l 



\w,0) 



(33) 



While the operator exp [i{(j)L ~ (pR)] corresponds to w = 1, and hence can have a non trivial one point function with 
N boundary conditions, the operator exp [i{(j)L — (j>R)/2] would correspond formally to w = 1/2. The corresponding 
state is not in the decomposition of the A'' boundary state, and therefore the one point function has to vanish. 

In a UV perturbation theory, the same sum over wells (which is compatible with the perturbation) will give rise 
to the same result. (Notice how the situation for the quantity exp [i(0H + '/>l)/2] is entirely different. In the UV, 
because the perturbation is cos $/2 = cos{(f)R + neutral combinations can be formed, and one finds non trivial 

pcrturbativc contributions. In the IR, the field obeys now (^/j + ^i,)/2 = (2n + l)n, and thus the exponential of 
interest is invariant when one goes from one well to the next.) 

Finally, the vanishing can also be understood in terms of form-factors. In the massive theory, the operator 
exp[i{(j)L — 4>r)/2] is not neutral, and increases the L ^ R charge by one unit. The non zero form-factors there- 
fore will have to involve an odd number of particles. Since the boundary state is a superposition of an even number 
of particles only, the one point function of this operator must vanish identically. A close correspondence with the 
disorder operator in the Ising model is discussed in the appendix. 

Of course, right at the UV fixed point, one can "by hand" restrict to one well of the potential, and find a non zero 
amplitude, of the expected form Gr{j)) <x It is tempting to force this amplitude to be non zero away from the 

fixed point as well. 

To do so, one can again use the form-factors approach. As explained in the appendix, the non zero value right at 

the UV fixed point can be obtained by adding to the boundary state a particle of vanishing energy (and momentum) . 
It is reasonable to try to define a non zero amplitude near the UV fixed point by generalizing the addition of this 
particle at finite value of the coupling. This leads to the result 



QR = {(ixp{i(j)L{-X2,y2)/2)exp{-i(j)R{xi,yi)/2)) 



1 



(.. - X. + ^iy2 - y^^^"^^'^ " + '^'^ " '^^^'^^'^ ^''^ 
where the second quantity has the form factor expansion 



i-(-)k - E ^ y ^ n I ^ tanh ^-e | J] (^tanh J 



(35) 



and the minus signs have been generated by the introduction of the particle at zero energy. 

Of course, this expansion is once again divergent. Moreover, in contrast with the case of the transmission amplitude, 
we do not even know the value in the IR, so we cannot divide by it to form a finite quantity. We can however observe 
that the logarithmic derivative of the correlator is regularizable by the same procedure previously used for {a)h- The 
first two terms read for instance 



d 
We 



\n{a{x,pB))h 



du 



PB 



-2xu 



-f 

Jo 



dudv , 



-(u -I- v)pb 



P'b 



(PB + u)^iPB +V)^ 



^-2x(u+v) _|_ 



(36) 



(where recall pB 



The integrals converge now at the origin, including when pb = 0- Call the series in the right 



hand side T,{xpb)- We can thus try to write 



1 r f^" dp' 

{a{x,pB))h = -y^exp / -^J:{xp'b) 

•"^ ' Uo Pb 



Pb 

Although right at = this expression reproduces the right result, there is a problem as soon as 7^ 0. This is 
because, at every order we have explored, T,{z) goes to a finite quantity as 2: — » 0, and thus the integral diverges at 
the origin! This can be studied in more details at the first order already, where the exponential reads, for small pb , 



exp 



p'b ^ 



To make the integral finite we have to introduce a new (IR) cut-off e, with dimension of temperature = length ^, so 
we get 



exp 



- ln(pB/e) 

TT 



(e/PB) 



1/tt 
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If we now let the cut-ofF go to zero, we see that the ampUtude again vanishes. A similar result was noticed independently 
by N. Sandler (private communication). 

This result is confirmed by higher order calculations. We think it is a consequence of the ill defined nature of the 
reflection amplitude within field theory, and can be interpreted by saying that this reflection amplitude is identically 
zero for any pB ^ 0, in the scaling limit. The only way to define this amplitude is to keep an additional (IR) cut-off, 
and then, to first order in the form-factors expansion, one finds 

It is in fact possible to obtain the exact value of the one point function by using more results about the Ising model, 
as detailed in the appendix: 

('^(^))/. = ^ M'^' (2*' - *) (38) 

where again 5* — ^'(l/2,l;z) = ^^Kq{z/2), z = 2pBX. From this the result follows that {a)j^ cx e^l'^x^l'^ at large 
values of psx. 

To each order we have explored, the value of the exponent at the origin was that is, every contribution but the 
first one seemed to vanish exactly. Meanwhile, the exact value expected from the formula H38() is 1/2. Most likely 
what happens is that the form-factors expansion does not converge uniformly at the origin - an identical feature takes 
place in the case of the ordinary one point function (cr)/i. 

Going back now to the physical amplitude we find 

Gi? cx -ij(ex)^/' VSP^IS*'-^') (39) 
In particular it follows that Gr{p) cx p^^^'^ (e/p)^^'^ in the IR hmit. 



dz 



z(l 



'2zxpB 



(37) 



D. Other amplitudes 

What one may call an "Andreev process" meanwhile would correspond to a Laughlin quasiparticle incident on the 
point contact giving rise to an electron in the same chiral Luttinger liquid, and a Laughlin quasihole in the other 
chiral Luttinger liquid. The amplitude for this process is 

^ ^g-i02(a;3>O,t3)/V2gi\/20i(x2>O,t2)g-i0i(3:i<O,ti)/V2^ ^^q^ 

The interacting part of this correlator, after using even odd combinations and folding, becomes 



and the argument has the form + 4>l — This time, the exponential is invariant when going from one well 

to the next, so the average will be non trivial. The amplitude vanishes both in the UV and IR limits, but has non 
vanishing perturbative contributions in the vicinity of either fixed point. All indications are that this amplitude is 
well defined in the scaling limit. 

In fact, the momentum dependent amplitude Ga{p,(i) (with say q the amplitude of the outgoing electron) has the 
same dimension as the would be amplitude for the reflection process (|39ll . It is likely that the regularized reflection 
amplitude we have defined previously H39|l coincides with the limit of Ga{p,<i) as q goes to zero (and then q has to 
be identified with e). 

Note finally that by folding one can map the geometry of two R moving Luttinger liquids to the geometry of a 
single non chiral Luttinger liquid. This can actually be done in several ways: the most natural is illustrated on the 
figure 12 In this case, what we called before transmission amplitude becomes a backscattering amplitude (small at 
small v) and what we called refiection amplitude becomes transmission amplitude (small at large v). 



III. POINT CONTACT BETWEEN A FERMI AND A LUTTINGER LIQUID 



A. Set-up 



A fundamental component in the experimental study of Luttinger liquids is the presence of a point contact with a 
three dimensional reservoir (in fact, a semi infinite three dimensional reservoir, whose boundary is a plane |3). The 
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TTie transmission am|:litude tiecomes a 
bad< scatering amplitude 



TTie reflection 
amplitude becomes a 
transmission amplitude 



FIG. 2: Tunneling between two identical chiral Luttinger liquids mapped onto the problem of a single impurity within a non 
chiral liquid. 



FIG. 3: Tunneling between two different chiral Luttinger liquids mapped onto the problem of tunneling between two different 
non chiral Luttinger liquids. 



analysis of the situation is somewhat similar to the Kondo problem. Electrons in the 3d reservoir can be organized 
not into plane waves but into spherical waves, and one finds that only the s — wave interacts with the Luttinger 
liquid (the Kondo impurity in the other case). Because there is a single mode, it can be considered as well as arising 
from a one dimensional conductor, and one can reduce the problem to the one of a point contact between a ID Fermi 
liquid and a ID Luttinger liquid. The situation then looks as in the figure |31 

By unfolding, one can map this problem onto the one of two chiral (say right moving) liquids with an interaction 
localized at the origin, corresponding to the bottom part of 13 The translation between the two is obvious (and 
different from the folding in figure |2Jl: for instance, a process where an electron comes in and bounces within the 
Fermi liquid corresponds in the unfolded picture to an electron going through on the top line. 

In a particular experiment, different tunneling terms will be induced by the realization of the point contact. At low 
energy - or low temperature - the decoupled fixed point, where the two electronic liquids are decoupled, is stable. To 
leading order, the tunneling that kicks in as energy is increased is a process where one electron from the Fermi liquid 
is exchanged for an "edge" electron in the Luttinger liquid. The amplitude for this term is traditionnally called T. In 
the limit of small T, the Lagrangian reads 



An 



dx)4>a + -^dx<j)b{dt - dx)<j)b 



rS{x) COs[-^(f>a - <j)b] 
\ V 



(42) 



where the field (/)a is associated with the Luttinger liquid, and (^^ with the Fermi liquid. This Lagrangian can be 
brought in a simpler form by the definition of new rotated fields ^a,h'- 



L = -^dxipa{dt - dx)ipa + -^dxfbidt ~ dx)(pb + T6{x) cos 
47r 47r 



1 



jiiPa - ^Pb) 



(43) 



where 



1 
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At high energy, tunneling tietween the two liquids is maximal 
> > 



At low energy, the two liquids are decoupled 

FIG. 4: Tunneling between the Fermi liquid and the = 1/3 Luttinger liquid as a function of the energy 



We will restrict here to = ^, which gives rise to g' — ^. The fields can now be gathered into an even combination 
which is free, and an odd combination that sees a boundary sine-Gordon type interaction. 

The perturbation induced by the tunneling term has dimension 2, and is irrelevant. This means that at low 
energies, there is essentially no tunneling between the initial Fermi and Luttinger liquids: one is at the fixed point 
where transmission is zero. This crucial feature of the problem, induces considerable complications. Indeed, suppose 
one were to consider an experimental realization starting from the trivial situation of decoupled wires at low energy, 
and increasing the energy to explore the non trivial scattering processed. In term of renormalization group, this would 
translate into going against the RG How, since the low energy fixed point is stable, and the coupling between the 
wires is an irrelevant perturbation. As is well known, such a situation does not have a well defined universal limit, 
as the RG trajectory is infiuenced by the many other irrelevant operators present in the problem. In our case, this 
means for instance that one also has to take into account the density density coupling, which also has dimension 2, 
as well as other fields maybe, and the results will be affected by the balance of these different terms. A priori, only 
extreme fine tuning would allow one to reach the perfect transmitting fixed point at high energy. 

Fortunately, there seems to be a way to realize the high energy fixed point T = oo. Indeed, in a very interesting 
paper, C. Kane pointed out that if tunneling between the two liquids takes place through an impurity, it is 
possible, by tuning two parameters only, to achieve perfect resonance that is, to have perfect tunneling between the 
Fermi and Luttinger liquids. Away from the resonance, deviations are controlled by the single relevant operator at 
the resonant fixed point, and the problem becomes formally identical to the by now well known problem of edge states 
tunneling at v = g' . 

In this new situation, at high energy, there is maximum tunneling between the Fermi and the Luttinger liquid which 
are strongly coupled, while at low energy, there is no remaining tunneling, the two are decoupled, and the scattering 
is entirely diagonal. This is illustrated on the figure^ 

The problem in the vicinity of the perfectly transmitting fixed point can be formulated in terms of dual fields, with 
Lagrangian 

L = -^dx^p+idt - dx)'f+ + -^dx^~{dt - dx)(p- + T6{x) cos ^/2g'(p- (44) 
47r 47r 

andf (X F-i/2. In this case, the perturbation is relevant, with dimension 1/2. As energy is lowered, physical processes 
where electrons and Laughlin quasiparticles are incident on the point contact and scatter into various combinations 
of outgoing electrons and Laughlin quasiparticles will have energy dependent amplitudes. 

In the following we denote the field (p^ simply by 0. We will also denote (pj^ by 0'. We fold the model so that 
instead of having right movers on the full line, we have left and right movers on the half line only, with the Lagrangian 

1 f° ~ 1 

^=^J dx[{dx^f -idt<^>f]+Tcos-^0) (45) 

In states are now right moving, and out states are left moving. As before, the perturbation has dimension 1/2, and 
is relevant. 

The most general quantities we are interested in read, in the original 4>a,b variables. 



G(4, 3, 2,1) — (T e^^P't>°a"*('^)gi'i'l>b''*(^)g~'--y^"'^T('^)g-^™'t'rW'^ 



(46) 
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m electrons q electrons 



n Laughlin q.p. p Laughlin q.p. 



FIG. 5: A process where m electrons and n Laughlin quasiparticles scatter into q electrons and p Laughlin quasiparticles. 



and are related (a precise definition would require, in particular, a choice of massless basis) to the amplitude for m 
electrons and n Laughlin quasiparticles to scatter into q electrons and p Laughlin quasiparticles, as represented on 
figure 13 

Following through the transformations detailed previously, this translates into 



G(4,3,2,l) = (T expi 



exp I 



(47) 



The fact that the (p' field is non interacting and sees Neumann boundary conditions imposes the (charge neutrality) 
constraint 

I + 9 = 3 + m (48) 

The full correlator will vanish in the limit of decoupled wires F = oo (F = 0) unlesss n — m = p — q and thus, 
q = m,p = n. This means that the scattering is diagonal, and that electrons and Laughlin quasiparticles go through 
without scattering, the well understood low energy fixed point. The correlator will vanish in the limit of totally 
coupled wires (perfect tunneling) F = (F = oo) unless n — m ^ —{p ~ q) . In this limit we have then, combining the 
two conservation equations, 

1 1 
q= -m + -n 

3 1 

p^-m--n (49) 

Note that at the high energy fixed point, strange things occur, as discussed in T|. For instance, all the processes 
where a single electron comes in m = l,n = 0, have vanishing amplitude! Still, probability must be conserved - so 
what may happen is that an infinity of processes have non vanishing amplitudes at non zero F, with the sum of the 
associated probabilities equal to unity, and that the corresponding series is not uniformly convergent: as F ^ 0, each 
individual term vanishes, but the sum remains equal to 1. We do not know how to explore this situation in more 
detail. 

For simplicity, we will mostly restrict to correlators of the form 

G(3,l) = g(3,l) (T exp{iV3a(l)'i^{-X3,t3))exp{--iV3a(j)'ji{xi,ti))) 

= 5(3,1)7 , ^. . TTT^, (50) 

where x^ > Q,xi < 0, and 

5(3,l) = (r ex:p{ia(f>L{-X3,t3))eTq){ia(t)ii{xi,ti))). (51) 

These correlators correspond either top — n — 0,q — 'moTq = m = 0,p = n and describe the scattering of q electrons 
into q electrons, or p Laughlin quasiparticles into p Laughlin quasiparticles. 

In the limit F = (high energy fixed point), the field ((> sees Neumann boundary conditions, and g vanishes by 
charge neutrality. This corresponds to the Fermi liquid and the Luttinger liquid being strongly coupled, so that all 
processes of finite amplitude involve non diagonal scattering, like electrons scattering into Laughlin quasiparticles. In 
the limit F — > 00, on the contrary, the field sees Dirichlet boundary conditions, the scattering of electrons or Laughlin 
quasiparticles is diagonal, and the correlator takes a finite value, 

5(3,1) -7 , _^ , ^ . . . , -TTT^' r^oo (52) 

In general, this correlator will exhibit a cross-over behaviour, which we will determine using the form-factors technique. 
Knowledge of the correlators will then allow us to determine the amplitude of interest. 
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B. The case 2e 



2e 



The case |a = 1| is special, and much simpler. Indeed, as discussed in the first part, the operator V± behaves like a 
fermion field, V± i oc tp^^"^ in a free fermion theory, the — soliton like a fundamental fermion, the + soliton like a hole, 



<0\^l\P>l- 

< 0\ijR\P >R- 



: /ie2 e 
/j,e 2 e 



-ie'^{t+x). 
ie'^{t-x). 



L+ < P\^l\0 >-- 



11 ^JU 



(53) 



where /i is a normalization factor (/i = v27r if one wants the fermions two point function to be normalized to unity). 
Similarly, for the dagger operators. 



It follows immediately that, for a = 1, using the same technique as in the first part. 



5(3, 1) =< 0|Vl(3)V1j(1)|S >=< V^L(3)^1^(1) >« 
Writing now the integral representation 

/•OO 

G(3,l)= / dp G(p,ps)e 
Jo 

we arrive, after convolution, at the following formula 



1 + e/3-/3B 27r 



p(xi-X34-i(j/l~J/3)) 



l! ^„e''(a;i-X3+i(yi-iy3)) 



such that 



G{p,pb) oc 



G{p,Pb) 
G{p,Pb) 



(p + Ps)^ln(l +p/pb) - 



p p 
6 ' 

PB^npB, 



3p^ + 2ppB 



(54) 



(55) 



(56) 



(57) 



(58) 



where we recall that the high energy limit is pb — >■ 0, F — > 0. 

Physically, the case a = 1 corresponds to q = m = 2, ie a process where a wave packet of charge 2e comes from 
the left in the Fermi liquid, and goes accross to a charge 2e wave packet still in the Fermi liquid. In the limit F = oo, 
this process occurs with probability 1, and it occurs with probability in the limit F = 0. 

The meaning of G as a probability in between these two limits requires some more discussion. The problem boils 
down to what one wishes to call a wave packet of charge 2e. In the bosonized theory, it is convenient to think of 
this as the object created by e^'^^-'^, which in the original Fermi liquid corresponds to the quantity x^^X (the fermion 
X is the fermion field in the Fermi liquid, and has no direct relation with the fermionic field which appears in the 
reformulation of the boundary sine-Gordon theory at the value — Att). When written in terms of modes, x^X reads 
as an superposition of states with two electrons at different momenta, weighed by some energy dependent coefficient. 
The corresponding creation operators have commutation relations of the form [a'' (k), a{k')] oc k^5(k — k'), and it is 
convenient to use these states as normalized states for charge 2e wave packets. With this convention, the probability 
that such a wave packet goes through unaltered is 



G[ip,PB) 
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with X 



G{ip, oo) 

This probability is represented in figure 



ix{l + ixY In 1 + - 



1 



3ix 
~2~ 



{ixf 



(59) 
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FIG. 6: Probabilities for the processes 2e — » 2e and e — » e. 



The case e 



The case a = 5, which corresponds to a process where an electron comes out as an electron, is of course much 
more complicated to study, and bears close resemblance to the transmission amplitude of the first part. A closed form 
expression is again possible, and one finds that, in imaginary time. 



G(3, 1) a /^ef-^/2^o(PBx/2), x = x^-xi^ t{y3 - yi) 

V 47r 



This provides a closed form expression for the Laplace transform 



G{p/pb) ^- f 

71" Jo 



dt 



which it is convenient to rewrite, to study the small pb limit, as 



1 rf 

G{p/pb) = -s/pbJp I 
Jo 



^t{l + t)^Jl-pBt/p 



(60) 



(61) 



(62) 



Prom this we see that the leading corrections goes in fact as y/pB/pln{pB/p) - in agreement with a perturbative 
calculation near the UV fixed point. The large pb limit can be studied as well, with 



(63) 



such that G = 1 when p/ps — 0. The leading order is linear in p/ps, again in agreement with perturbative calculations 
near the IR fixed point. It is convenient to reexpress G as a hypergeomctric function 



G(p/pb)=F(1 i 1;- 

\2 2 Pb J 



P_\ 



(64) 



Standard formulas then allow for a complete expansion in the large pb limit. For further use, we quote the general 
result 



W2)^ai,c;-^ 



r(c) 



(Pb/p)^/' gr(l/2 + n)r(3/2-c + n) 1 



2'2'' pbJ r{c-l/2)^^ r(l/2) r(3/2-c) (n!)H fB 



P 



[\n(p/pB) + hr,{c)] (65) 
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where 

K {c) = 2?A(1 + n)~ V(" + 1/2) - ?A(c - n - 1/2) (66) 

In the present case, G, after analytic continuation, can directly be interpreted as the scattering amplitude for an 
incident electron into an outgoing electron, and |Gp is therefore the probability of this process 



Pl_l = \G{ip/pB)f = 



,11 n ^ ^ 



(67) 



The result is also shown in figure El 

Note that we can also extract from these results expressions for the density of states inside the Fermi liquid. This 
density has a constant part and an oscillatory part, and thus reads, as a function of the distance x from the point 
contact: 



p(a;) cx est + e*P^^ y dp G{ip/pB)e'P'' + h.c. (68) 
This will be discussed in more details elsewhere. 

D. The case Iqp Iqp 

When a Landau quasi-particle (Iqp) scatters into a Landau quasi-particle, the amplitude is related to a slightly 
different kind of correlator, 

G(4,2) = 5(4,2) (T exp(^(a/y3)0'i(-X4, ^4)) exp(-z(a/V3)(/.i^(a;2, ia))) 

= ^^(^'2)7 , ^, . f, TTwTI' (69) 

(X4 — X2 — 14 + t2 + le sign{t4 — 12) r ' 

where 2:4 > 0, 2:2 < and 

5(4,2) = (T exp(-iQ;(/)L(-a;4,t4))exp(-iQ;0_R,(x2,t2)))- (70) 

The case Iqp — > Iqp corresponds again to q = 1/2, that is an expression which is identical, in the (j) sector, to the 
amplitude for the electron electron process. One can then deduce the Ipq amplitude from the electron amplitude 

roc /•! / \ —1/2 

G(4,2)cx / p-^/^e'f"=dp {1 - uy^/^u-^/^ ( 1 + ^u] du (71) 
Jo Jo \ PB J 

where a; = 2:4 — a;2 + i(jj4 — 1/2)- To obtain the latter result, we used a Laplace representation of x^^^, which is 
divergent, and can be made finite by the analytical continuation (continuation of a F function in this case). Similarly, 
(|71|l is divergent, and can be made finite by expanding in powers oi p/ps at low energy, and using B ( F) functions 
rcgularization. One finds then 

t^o^P^^ 22"(n!)2 F(n+l/3) 

By applying Stirling's formula at large n, one can check that the series converges for p/ps < 1- At large energy 
meanwhile, one finds the leading behaviour 

-2/3 I Pb} 



G{p,PB)^p-'"\^^j Hpb/p) (73) 

again in agreement with perturbative calculations. The complete large energy expansion can be obtained by recog- 
nizing that 



(74) 



G(„a)«.-"f(l.l,i;-A) 

and then using (|^ . 

A graph of the associated probabihty (normalized by the low energy limit) i.e. \F{^, ^, provided in 

figure (O. 
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E. The case le 



This is still another case, corresponding to m = l,p = 3. The amplitude is therefore 

1 



G(4,l) = 5(4,1) 



where X4 > 0,xi < and 



{x4 — Xi — + ti + ie sign{t4 — ti))"^/'' ' 



3i 



5(4,1) = (T exp -—(l)Li-X4,t4) cxp -(f)R{xi,ti) ]). 



(75) 



(76) 



This amplitude is expected to vanish both at the UV and IR fixed points, and take non trivial values in between. Note 
that the amplitude is invariant under the shifts (pR ^ (pR + w and (f)L ^ (f>L — corresponding to the UV boundary 
state, and 4>r ^ (t>R + 27r, (f>L ^ (f>L + Stt corresponding to the IR boundary state. One can similarly check that the 
leading relevant (resp. irrelevant) operator give it a non trivial value away from the UV (resp. IR) fixed point. More 
sophisticated arguments presented in the appendix give rise to the following expression 



G(4, 1) oc x-^ r e-y{l - 2y)y-^'^ f 1 + ^ 
Jo \ Pbx 



-1/2 



dy 



from which the momentum dependent amplitude follows 



G{p,pb)(xp 



1/4 



[\i-ur 

Jo 



1/4^-1/2: 



1-w 



1 + —u 1 du 
Pb 



P 



(77) 



(78) 



Note that this correlator vanishes in the low and high energy limits, as expected from previous considerations about 
the fixed points. An expansion at low energy follows from T function regularization again 



G{p,pB)o,p'/^^{-iy 



(2n)! 



22"(n!)2 



r(-l/4)r(l/2 + n) r(-l/4jr(3/2 + n) 



The leading behaviour at high energy is 



r(l/4 + n) 



r(5/4 + n) 



P_ 
Pb 



(79) 



(80) 
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FIG. 8: Probability (arbitrarily normalized) for the process e — » Cedge- 



Once again, one can recognize a sum of two gamma functions 

"r(i/2)r(3/4) „ /I 1 5 



r(5/4) pb) 



r(i/2)r(-i/4) 
r(i/4) 



F 



1 1 1 

2' 2' 4' 



PB 



(81) 



and then use formula H65fl to obtain the complete high energy expansion. 

In the present case, we do not know how to normalize the associated probability, which is expected to vanish in 
both low and high energy limits - this would presumably require a more complete understanding of unitarity issues. 
Qualitatively, we expect that it should behave as \G{ip,pB)/p^^'^\'^, which is represented in figure 



IV. CONCLUSIONS 



The results we have obtained in closed form are but a very small step towards a general understanding of the 
scattering amplitudes. Extensions of this work to other problems would involve in general considerable technical 
difficulties, if only because the form-factors approach is so hard to implement in the massless case. On top of this, 
even the form-factors themshelves might not be known, or only known in a very involved form - clearly, a new idea 
is needed there if progress is to be made. 

A more rewarding direction of development would probably be to restrict to the present g = \ case, and consider 
instead the question of unitarity, and how the probability spreads among multiparticle processes as energy is varied. 
It also remains to be seen whether a direct determination of scattering amplitudes is experimentally feasible. 
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APPENDIX A: DISORDER OPERATOR AND FLOW CONSIDERATIONS. 



The one point function of the disorder operator in the Ising model with a boundary magnetic field vanishes all 
along the flow. This follows immediately from the fact that all non vanishing form-factors for ^ involve odd numbers 
of particles, while the boundary state is a superposition of an even number of particles only. 
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Here, a word of caution is necessary. Right at h = 0, the proper boundary state in the massive theory is in fact a 
superposition of states with odd and states with even numbers of particles: 



\B> 



free' 



[l + Z\9)) exp 



Kfr,Mz^{-e)zHe) 



(Al) 



One has Kjree = — icoth|, and the presence of Z^O) is related with the pole of K at 9 — 0. So, right at the free 
boundary conditions fixed point, the one point function of the disorder operator is non zero. In the massless limit, 
the zero momentum particle becomes a zero momentum zero energy particle, and becomes almost invisible except for 
the fact that it changes odd into even number of particles, and intorduces some sign factors. One then finds (this is 
discussed in more details below) 



n!./__AA 1 2n 



Yl ( tanh 



(A2) 



where the term Yli tanh which is there equal to (—1)" , is cancelled by the (—1)" term arising from the presence 

of an additional particle at vanishing energy and momentum. This formal expression coincides with {(y{x))h=oo, and 
thus {ii{x))h=o oc x~^^^. 

For any non zero h however, there is no one particle state in |i? > jll^ . and the one point function of /i vanishes 
exactly. 

In this appendix we first would like to discuss this result further. To do so, we first can try to follow the route of 
0, and derive a differential equation for < fj, >. Recall that in the usual case, Chatterjee and Zamolodchikov find 
the following equation 



d 

dY 



1 

2Y 



16Y^ 



(cr) =0, r = 2pBX 



{x taken positive by convention here) and pb = 'iirh^. Setting (cr) = y-'/^vj/^ one finds the equation 

Y^" + (1 - Y)^' - ^ = 



(A3) 



(A4) 



whose physica Isolution is proportional to the degenerate hypergeometric function — ^'(1/2, 1,1"). 

In the case of the disorder operator, the same logic as in [g can be followed, with proper modifications due to the 
different monodromies involves. We find instead the equation 



r d^ 




d 


(-^ 


['dY^ + 




dY + 


V 2Y 



1 



16^2 



(pl) = 0, r = 2pBX 



And, setting (fi) = Y ^Z^^*, one has this time 



y^r" _ Y^' = 



(A5) 



(A6) 



with only solution (that does not grow exponentially at large distance) ^ a constant. The only physical solution is 
thus < /i 0, since we know it has to vanish in the limit of large fields, where the boundary spins are fixed (the 
case pb — can be exceptional, since then the variable in the differential equation is not well defined). 

This calculation is formally the same as the calculation of the one point function of the spin (cr)^ for an Ising 
model where the boundary perturbation, instead of being the boundary field, would induce the following equations 
of motion: (if t is the coordinate along the boundary as in p3]) 



(A7) 



Note the crucial minus sign on the right hand side - in the usual problem of Ising model with a boundary magnetic 
field, the rhs comes with a plus sign, and the physics is the flow from free to fixed boundary conditions. It is easy to 
check that switching the sign in this equation leads to switching the sign of the reflection matrix, which reads in this 
case 



R 



-i tanh 



in 
T 



(AS) 



(in the following we set now pB = ATrfi^). In the UV, one finds R ~ ~i and in the IR i? = i, so the role of the free 
and fixed boundary conditions are exchanged: this is the flow 'dual' to the ordinary one, which can be understood 
for instance by using a low and high temperature graphical expansion of the Ising model partition functions. 
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Here, one has to be a bit careful. From entropy considerations, one does not expect a flow from fixed to free boundary 
conditions to be possible. Rather, what probably happens is a flow from the superposition of two boundary states 
(fixed + and fixed — ) towards free boundary conditions. The ratio of degeneracies is then guv/gm ^ ^ x = V2, 
the same as for the usual flow from free to fixed. In such a flow, the Z2 symmetry is never broken, and therefore 
< (T >^= always. The free energy in this flow is in fact the same (as a function of T,pb) as the free energy in 
the ordinary flow - a feature compatible with the fact that the free energy in the integrable approach only depends 
on the derivative of the log of the R matrix. In this point of view the UV fixed point appears extremely singular, 
as it is a superposition of two independent states - for any ps 7^ 0, these two states are coupled, and "cutting the 
system in {Z2) half" does not make physical sense. Only at the "double fixed" fixed point can one do so, and define 
(^>^=o « a;-i/8. 

Now we can ask about the boundary state associated with (|A8p . In the usual case, the proper way to discuss this 
is to go to the massive theory. This can be done here in various ways. An interesting one is to make a detour through 
the massive Kondo model |24j . that is the theory of massive fermions in the bulk, with a Kondo like interaction at 
the boundary. This theory decouples into two independent Ising models. One of them is the usual Ising model with 
a boundary magnetic field, while the other has the following R matrix 



ITT u \ 1 + + i sinh 9 



hi 

R^-i coth I — - 1 1 ' • ' (A9) 

2m 

Here one sees that R — Rfree at large value of the coupling h, while R — R fixed at small value of the coupling. For 
any value of the coupling, the K matrix has a pole at the origin, indicating the presence of a zero momentum particle 

in the boundary state. Note that in the massless limit the two reflection matrices become R ~ i tanh ( 



and R — —i tanh — . The model thus decomposes into two Ising models, one with a flow from free to fixed, 

and one with a fiow from "double fixed" to free (note that the ratio of g factors is thus guv / giR = \/2 x = 2 as 
required for the spin 1/2 Kondo model.) The boundary state for the latter flow always contain a zero momentum, 
zero energy particle 

To proceed, let us consider the form factor 

2n 2n 

(O|^|O,02n,...,^i) =ntanh(^./2) H tanh(%/2) (AlO) 

i—l i<j — l 

In the massless limit, n particles will become R movers, n will become L movers. Whatever particles we choose for 
this, the particle at rapidity zero will introduce a factor (—1)", and multiplying this factor, we will have the massless 
limit of the product nj<7=i tanh(0y /2), which will be a factor of similar terms for R and L components independently, 
times a sign factor which depends on which particles become R and which particles become L. The minus sign hence 
generated by the particle at vanishing rapidity exactly cancels the minus sign coming from the R matrix, so we can 
write at the end 

= (AH) 

hence establishing that the one point function of the disorder operator in the dual flow is the same as the one point 
function of the order operator in the original fiow. 

A last comment about the dual fiow. Although we do not know the action corresponding to the one point function 
near the fixed boundary conditions fixed point (and thus the parameter h is only formal), we can find it near the IR 
(free) boundary conditions fixed point by using the arguments of 0- Since the R matrix is analytic, it follows that 
this action involves only terms of the form ^pd'^^'^tp. Since on the other hand the correlation functions of the spin 
operator with these (conserved) quantities all vanish with free boundary conditions, it follows that, in a perturbative 
expansion near the free boundary conditions (IR) fixed point, (cr)^ has to vanish identically. Of course, it may well 
be that the expansion of this quantity has a vanishing radius of convergence, wiht non perturbative contributions - 
the foregoing arguments have established that this is not the case however. 

As discussed in the main text, a regularized version of {ij)h can be obtained by adding to the boundary state a 
fermion at vanishing energy and momentum. The corresponding expression in terms of form-factors is given in Ij^SI) . 

An exact expression for this can be obtained by considering, following Chatterjee Zamolodchikov j^, the correlator, 
< x(2:)/i(w, w) > where x = (dz + iXjip for Im z > which we assume in the following. To start, we consider in fact 
a simpler correlator, < 'il){z)fi(w,w) >. We calculate this correlator using the form factor formalism, after insertion 
of the boundary \B >. The n*'' order term then reads 

^ / n df3^{-i) tanh ^^^(-l)"("-i'/2 ^ o|^(z)/i(u;, u;)|/3i, ...,/?„; /?„ >r.,...,r,l,. 
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where the (— originates from reordering the n*^ temr in the expansion of the boundary state. Note that 
this term is exactly cancelled by the sign coming from the massless limit of the ^ form factor. 

The disorder operator has non vanishing form factors only on odd numbers of particles. The fermion is simply able 
to destroy one particle. Suppose /i destroy the particles (32, ■ ■ ■ , in the R channel and . . . , /3„ in the R channel. 
Then ip the remaining particle f3i in the R channel, so the matrix elements < 0|... > will give the contribution 



=/3i/2 



n 

Y[ tanh 

i<j=2 



JJ^ tanh 



(3k - (3i 



exp [—2a; (e 



)] exp[-(z- w)e^i] 



k<l=l 



and X = 2£±H. 



Let us now look at the large z behaviour (/m z > 0) of the corresponding integral. In that limit, the integral is 
dominated by the region where /3i — oo, and we are left with the integral over P2, ■ ■ ■ Pn of the expression 



^" ^ tanh^?^^(-l)' 



J] tanh ^i-^ exp [-2a; (e^^^ 



0] 



i<j=2 



k<l=2 



while the integral over /3i produces a factor \l \pz. The crucial factor (—1)" comes from the limit of the form factors 
contributions tanh . The factor coming from the boundary state meanwhile reduces to n"=2 tanh . 

Now we could also have "0 annihilate any of the other /32, .../?« R particles instead. If this so happens, commuting 
-0 through the particles on the left of the ones it annihilates generates minus signs. These however will be cancelled 
exactly by the minus signs coming from the limit of the product of tanh factors. It thus simply follows that we get a 
multiplicity of n for the term we have just analyzed, and therefore at large z 



00 /-oo ^ dj^- / 



pB-(3^ 



(-1)" Yl (tan^ 



A - A- 



exp 



~2x {e^^ + . . . + e"-)] 



and the proportionality factor is independent of A. Of course as always the generic terms in the series are divergent, 
and we do not learn much (in particular, the z power law dependence is not reliable) until we regularize by taking 
ratios. 

Using , we can on the other hand calculate the large z behaviour by computing B in the right hand side of their 
equation (24). Putting things together we obtain the final result that 



n=0 •^-°°i=l 



Y[ (tanh 

i<j=i 



^3^ - 13, 



exp [—2a; I 



)] 



/ (same with 13b = 00) = — \/z 



2— vI/(l/2,l;z)-vE'(l/2,l;z) 
az 



2pBX 



(A12) 



Numerical checks of this result (obtained by carrying out the form factors expansion to fourth order) are presented 
in the figure. The agreement is not spectacular, but maybe this is expected, as successive orders in the form-factors 
expansion all add up with the same (positive) sign, while in the usual case they form an alternating series. Note the 
divergence at vanishing pB- 



APPENDIX B: ANALYTICAL RELATIONS. 

This appendix details a simple extension of the method of Chatterjee-Zamolodchikov to determine correlators 
directly in the sine-Gordon model. Let us start from 

1 r° 1 

dx[{d^^f -{dt^f\+vcos-^{x = Q,t). (Bl) 

Fermionizing using (we revert to the more standard notation V':V')'0±(^) = e^*'^^^-' and 'ji+jz] = e^^'^^^K we get an 
action consisting of two decoupled Ising models having an interaction at the boundary (see |23 for details). Taking 
the variation with respect to the boundary action, we find boundary conditions for the fermionic fields 

0+(z) + 0_(z) = 0+(z) + 0_(z) (B2) 
idt{ip+-ip-) + v{ip+-'ip-) = -idt{'ip+-il^-) + v{Tp+-7p_). (B3) 
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FIG. 9: Exact result and form factors estimates for the Ratio ljA12|l as a function of ps. 



Following Chatterjee-Zamolodchikov, we define then 



(B4) 
(B5) 



with z = ^{t + ix). The boundary conditions insure that the fields are analytic in the full complex plane. In terms 
of the bosonic field they read 



X+{z) = 2cos(/)(z) 

X-{z) = — 9z sin 0(z) + if sin 

Now, given the OPE of the operators x± with the vertex operators V\/2,-i/2 we expect that 



{z — w) 



(B6) 
(B7) 

(B8) 
(B9) 



(z — w) 

The first equation for instance follows from analyticity together with the OPE around z = w 

X+{z)Vi/2^^i/2{w,w) ~ ^^_^^y/2 (BIO) 

+ iz~w){V3/2,-i/2-idW-i/2.-i/2 + ---)] (Bll) 

a similar OPE around z = w, and the condition x+ ~ s-t infinity. 

Expanding the first equation in (|B10|) in powers of z — w leads to the determination of -6+ = (w — 
u))^/^(y_i/2,-i/2(w, iS)). The comparison with the expansion in (z — w) leads to a tautology relating the one point 
fimction of 14/2,1/2 ^-nd that of l^_i/2,-i/2- 

Expanding to higher orders gives non trivial relations. For example, to next order in (z — w) we have 

B{w, w) 



2{w - w)3/2 



(^3/2,-l/2> - i(50^-l/2 -1/2) 



leading to the simple 



1 d 

(V^i/2,-l/2> = (VLl/2,-1/2) - 2— (F_i/2,_i/2) 



(B12) 



(B13) 



giving an exact result for {V^^2,-i/2) in terms of (VLi/2,-i/2)- Since the latter quantity is known from |8(|, the former 
follows. 
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